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Summary. 

The  problem  of  a  semi— infinite  mode  III  crack  that  suddenly  begins  to 
propagate  at  a  constant  speed  is  considered  for  a  general  linear  viscoelastic  body. 
Simple  closed  form  expressions  are  derived  for  the  Laplace  transform  of  the  stress, 
displacement,  and  stress  intensity  factor.  A  Barenblatt  type  failure  zone  is  assumed 
to  exist  at  the  crack  tip  and  an  expression  for  the  total  energy  flux  into  the  failure 
zone  per  unit  crack  advance,  G(t),  is  derived.  The  short  time  asymptotic  form  of  G 
is  constructed  and  the  rate  at  which  it  converges  to  its  steady-state  limit  is  studied. 
It  is  shown  that  the  rate  of  convergence  to  steady— state  is  dependent  upon  crack 
speed  and  material  properties.  Moreover,  it  is  found  that  whether  or  not  a  failure 
zone  is  incorporated  into  the  model  significantly  influences  both  quantitatively  and 
qualitatively  the  short  and  long  time  behavior  of  G.  This  difference  is  important  to 
predictions  of  stable  vs  unstable  crack  speeds  based  upon  a  critical  work  input  to 


the  failure  zone  fracture  criterion. 
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1.  Introduction. 

Several  analytical  studies  of  dynamically  propagating  cracks  in  linear 
viscoelastic  material  have  appeared  in  the  literature  since  Willis(  1967 )  presented  an 
analysis  of  the  dynamic,  steady— state  propagation  of  a  semi— infinite,  mode  111 
(anti-plane  shear)  crack  in  an  infinite  viscoelastic  body.  Employing  transform 
methods  and  the  Wiener— Hopf  technique,  Willis  constructed  the  dynamic  stress 
intensity  factor  (SIF)  for  a  standard  linear  solid  material  model  and  general  crack 
face  loadings.  Subsequently,  Atkinson  and  List(1972)  introduced  transient  effects 
into  the  problem  by  assuming  that  the  crack,  initially  at  rest,  begins  to  propagate  at 
a  constant  speed  under  the  action  of  suddenly  applied  loads  on  the  crack  faces.  Also 
utilizing  the  Wiener— Hopf  method,  they  derived  an  expression  for  the  Laplace 
transform  of  the  time  dependent  SIF  from  which  long  and  short  time  asymptotic 
approximations  and  numerical  Laplace  inversion  calculations  were  obtained. 
However,  their  analysis  was  limited  to  consideration  of  constant  applied  load  along 
the  crack  faces  and  the  required  Wiener— Hopf  factorization  was  effected  only  for  a 
Maxwell  material  model  and  the  Achenbach— Chao(  1962)  three  parameter 
approximation  to  the  standard  linear  solid. 

Somewhat  later,  Atkinson  and  Coleman(1977)  used  a  matched  asymptotic 
expansion  technique  to  develop  an  approximate  analysis  of  the  steady— state 
propagation  of  a  semi— infinite  mode  I  (plane  strain  or  plane  stress)  crack 
propagating  in  a  clamped  viscoelastic  strip.  Shortly  thereafter,  Atkinson(1979) 
presented  an  approximate  analysis  of  the  mode  I  counterpart  to  the  mode  III 
problem  considered  by  Atkinson  and  List(1972).  Their  argument,  ostensibly  valid 
for  fairly  general  material  models,  involved  slightly  modifying  the  exact  elastic 
result  of  Baker(1962)  in  order  to  approximate  the  Laplace  transform  of  the  actual 


dynamic  viscoelastic  SIF.  The  dominant  term  for  each  of  the  short  and  long  time 
asymptotic  expansions  of  this  approximate  SIF  was  derived  for  each  of  three 
different  applied  crack  face  loads:  a  constant,  a  delta  function  point  force  and  an 
exponentially  decreasing  form.  Also  in  that  paper,  Atkinson  reconsidered  the  mode 
III  problem  and  extended  the  Atkinson  and  List  analysis  to  handle  the  above  three 
types  of  crack  face  loadings.  However,  consideration  was  limited  again  to  the 
Achen bach— Chao  material  model.  Atkinson  also  constructed  an  expression  for  the 
energy  release  rate  (ERR)  based  upon  a  local  (i.e.  at  the  crack  tip)  work  argument 
and  the  singular  stress  field. 

Also  in  that  same  year,  Atkinson  and  Popelar(1979)  presented  an  analysis  of 
the  transient  constant  crack  speed  mode  III  problem  for  a  viscoelastic  strip. 
Constitutive  relations  in  terms  of  differential  operators  were  assumed  and  the 
external  load  consisted  either  of  constant  displacement  of  the  upper  and  lower  layer 
boundaries  or  constant  tractions  on  them.  The  crack  faces  were  assumed  to  be 
stress  free.  Again  the  Wiener— Hopf  method  was  used  to  construct  an  exact 
expression  for  the  Laplace  transform  of  the  SIF.  The  required  Wiener— Hopf 
factorization  was  carried  out  modulo  a  term  involving  a  Cauchy  type  integral. 
Atkinson  and  Popelar  then  restricted  attention  to  numerically  approximating  the 
Cauchy  integral  for  the  steady-state  limit  case  and  assumed  a  standard  linear  solid 
material  model. 

A  year  later,  Atkinson  and  Popelar(1980)  addressed  the  more  difficult  mode  I 
problem  for  a  viscoelastic  strip.  Again  by  use  of  the  Wiener— Hopf  method,  a  formal 
expression  for  the  Laplace  transform  of  the  SIF  was  constructed  containing  a 
complicated  Cauchy  integral.  As  with  the  corresponding  mode  III  case,  the  integral 
was  studied  numerically  in  the  limiting  special  case  of  steady— state  crack 


propagation  in  a  standard  linear  solid. 

Somewhat  later,  Walton(198'2)  examined  further  the  steady-state  mode  III 
problem  considered  by  Willis(1967).  Utilizing  the  Riemann— Hilbert  rather  than  the 
Wiener— Hopf  methodology,  he  constructed  a  simple  closed  form  expression  for  the 
SIF  valid  for  general  crack  face  loadings  and  very  general  material  models.  More 
specifically,  constitutive  equations  expressed  in  terms  of  convolution  integrals  rather 
than  differential  operators  were  adopted  and  the  results  were  shown  to  be  valid 
irrespective  of  any  assumed  time  rate  of  decay  of  the  viscoelastic  shear  modulus.  In 
contrast,  constitutive  relations  in  terms  of  constant  coefficient  differential  operators, 
necessarily  force  an  exponentially  decaying  modulus  thereby  preventing 
consideration  of  the  important  class  of  power— law  models  which  more  effectively 
represent  the  mechanical  response  of  many  real  viscoelastic  materials,  such  as 
rubber,  than  do  exponentially  decaying  functions. 

Subsequently,  WaIton(1983)  extended  the  above  analysis  to  determine  the 
angular  dependence  of  the  stress  field  in  a  neighborhood  of  the  crack  tip.  In 
particular,  it  was  shown  that  the  asymptotic  stress  field  at  the  crack  tip  has  the 
same  angular  dependence  as  the  corresponding  dynamic  elastic  problem.  Only  the 
SIF  differs  between  the  elastic  and  viscoelastic  fields. 

Walton(1985)  next  considered  the  steady-state  mode  III  problem  for  a 
viscoelastic  strip.  Again  utilizing  the  Riemann— Hilbert  method,  a  closed  form 
expression  for  the  SIF  was  constructed  for  general  loadings  and  shear  modulus.  The 
form  of  the  solution  exhibits  clearly  the  combined  effects  of  material  properties, 
crack  speed  and  layer  thickness  upon  the  SIF. 

More  recently,  Walton(  1987a)  reconsidered  steady-state  mode  III  crack 
propagation  in  an  infinite  viscoelastic  body  in  order  to  investigate  the  implications 
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of  including  a  failure  zone  model  of  Barenblatt(1962)  type  into  the  determination  of 
the  ERR  from  a  global  energy  balance  calculation.  A  simple  closed  form  expression 
for  the  ERR,  which  in  this  case  is  just  the  work  done  by  the  tractions  in  the  failure 
zone,  was  derived  under  the  same  mild  conditions  on  the  shear  modulus  assumed  in 
Walton(1982)  and  for  a  fairly  broad  class  of  crack  face  and  failure  zone  loadings.  It 
was  then  observed  that  whether  or  not  a  failure  zone  is  incorporated  into  the  model 
greatly  influences  both  qualitatively  and  quantitatively  the  dependence  of  the  ERR 
upon  crack  speed  and  material  properties.  In  particular,  calculations  based  upon  a 
failure  zone  seem  to  reflect  more  closely  experimental  observations  of  cracks  rapidly 
propagating  in  real  viscoelastic  material.' 

The  methods  of  Waltonf  1987a)  can  be  applied  to  the  calculation  of  the  ERR 
for  a  wide  variety  of  dynamic  viscoelastic  crack  problems.  Schovanec  and 
Walton(  1987c)  recently  completed  the  analysis  of  the  dynamic  steady-state 
propagation  of  two  parallel  mode  III  cracks  in  an  infinite  viscoelastic  body.  Also, 
Walton(1987b)  has  recently  completed  a  study  of  the  mode  I  analog  of 
Walton(1987a).  In  both  of  these  investigations  a  Barenblatt  failure  zone  model  was 
adopted.  Of  related  interest  are  two  additional  papers  by  Schovanec  and 
Walton(  1987a, 1987b)  in  which  these  same  methods  were  applied  to  quasi-static 
mode  I  crack  propagation  in  non— homogeneous  viscoelastic  material.  It  should  also 
be  noted  that  Knauss(1973)  and  Schapery(1975)  applied  the  Barenblatt  model  to 
quasi-static  viscoelastic  crack  growth  in  viscoelastic  material  and  observed  that 
whether  or  not  a  failure  zone  is  incorporated  into  the  model  greatly  affects  the 
behavior  of  the  ERR. 


I  -  — - 
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1  Private  communication  with  Prof.  W.G.  Knauss  and  Prof.  K.  Kuo. 
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The  present  paper  applies  the  above  Barenblatt.  failure  zone/ 

Riemann— Hilbert  program  to  the  transient  mode  III  problem  of  a  semi— infinite 
crack  in  an  infinite,  linear  viscoelastic  body  that  begins  to  propagate  at  a  constant 
speed,  v,  due  to  the  sudden  application  of  tractions  to  the  crack  faces  which  then 
propagate  with  the  moving  crack  tip.  Thus  the  crack  speed  is  constant  for  times 
t>0  and  the  crack  face  tractions  are  spatially  constant,  though  time  varying, 
relative  to  the  moving  crack  tip.  Clearly  it  would  be  desirable  to  allow  nonsteady 
crack  growth  in  order  to  study  the  initial  crack  acceleration  phase.  However,  for 
dynamic  crack  propagation  in  viscoelastic  material  such  an  analysis  seems  quite 
difficult  to  carry  out.  Moreover,  the  primary  issue  addressed  here  is  the 
determination  of  the  influence  that  crack  speed  and  viscoelastic  constitutive 
properties  have  upon  the  rate  of  convergence  to  steady  state  of  the  work  input  to 
the  failure  zone.  Consequently,  it  is  felt  that  valuable  insight  into  the  asymptotic 
behavior  of  the  more  realistic  model  with  unsteady  crack  growth  can  be  obtained  by 
studying  the  admittedly  artificial  transient  problem  in  which  the  crack  speed  is  kept 
constant.  For  the  sake  of  mathematical  completeness,  a  short  time  asymptotic 
analysis  is  also  presented  here. 

The  next  section  contains  a  derivation  of  the  asymptotic  singular  stress  field 
in  front  of  the  advancing  crack  tip  in  the  absence  of  a  Barenblatt  failure  zone.  More 
specifically,  a  closed  form  expression  for  the  Laplace  transform  of  the  stress 
intensity  factor  is  constructed.  The  analysis  is  valid  for  a  very  general  class  of  shear 
moduli,  p( t),  and  a  very  general  class  of  time  and  spatially  varying  crack  face 
tractions,  thereby  generalizing  the  results  of  Atkinson  and  List(1972). 

In  section  3,  a  Barenblatt  failure  zone  is  introduced  and  an  expression  for  the 
work  input  to  the  failure  zone  is  derived  for  general  shear  modulus  pi t)  but  for 
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special  crack  face  applied  tractions  and  failure  zone  constitutive  model.  The  failure 
zone  constitutive  model  adopted  here  is  certainly  open  to  criticism  on  physical 
grounds.  But  again,  for  reasons  discussed  in  section  3,  it  is  felt  that  consideration  of 
this  admittedly  artificial  model  provides  valuable  insight  into  the  longtime 
asymptotic  behavior  of  more  realistic  models  under  the  assumption  of  small  scale 
yielding.  The  final  section  contains  the  long  and  short  time  asymptotic  analyses. 
The  results  are  then  illustrated  through  special  cases  and  examples. 
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2.  Problem  Formulation  and  Stress  Analysis. 

The  problem  to  be  considered  is  that  of  a  semi— infinite  mode  111  crack  that 
begins  to  propagate  at  a  constant  speed  v  in  an  infinite  viscoelastic  body  due  to  the 
sudden  application  of  crack  face  tractions  that  then  travel  with  the  crack.  The 
governing  field  equations  for  the  motion  of  a  linear  viscoelastic  solid  are 
pu.  =  £r..  -+U'  •  ),  <r..=2u*de. .  +  &.A*de,  ,  , 

where  <7p  e-,  and  u-  denote  the  stress,  strain,  and  displacement  fields  respectively. 

In  (2.1),  fi*d(  denotes  the  Riemann— Stieltjes  convolution  /i*de=  f4t— r)  de(r). 

Since  the  deformation  is  assumed  to  be  antiplane  strain,  ui=0,  U2=0,  and  the 
only  equation  of  motion  not  identically  satisfied  is  /i*dAu3=pu3,  where  A  denotes 
the  Laplacian  operator. 

A  semi-infinite  crack  lying  along  the  negative  X|-axis  is  assumed  to  begin  to 
propagate  at  time  t=0  with  a  constant  speed  v  driven  by  loads 
^23(xi,0,t)=L  A  (x-~vt,t)  which  follow  it,  where  A  (  )  is  dimensionless  while  a  and 

6  6  3,  6  c 

e 

Lg  have  the  dimensions  of  length  and  stress,  respectively.  Thus,  while  the  crack 
speed  is  assumed  to  be  constant,  the  driving  load  is  allowed  to  be  time  varying.  The 
corresponding  initial— boundary  value  problem  is 

/>U3=^*dAu3=A|p(0)u3(xi,X2,t)+J^  us(xi,X2,t)  p'(t-r)  drj 

with  initial  conditions  u3=0,  u3=0  at  t=0 

and  boundary  conditions  ^23(xi,0,t)s=L  A  ,t)  X(<vt 


e  ev  a. 


u3(xi,0,t)=0 


X|>Vt 


(2.1) 


(2.2) 


(2.3) 


2  2 

^23(xj,X2,t)-*0  as  Xi  +X2  -*0O. 


From  (2.1)  it  can  be  seen  that 
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<?23(X|,X2,t)=|^2(/t*dU3)=^jK0)u3(X|,*2,t)+J^U3(x,,X2,T)/i'(t-T)  drj. 

It  is  convenient  to  change  from  the  fixed  coordinates  (xi,x2,t)  to  the  moving 
coordinate  system  (x,y,t)  given  by  x=xj— vt,  y=x2  and  to  define  w(x,y,t)  by 
w(x,y,t)=w(xi— vt,X2,t)=U3(xi,X2,t).  In  the  moving  coordinates,  equation  (2.2) 
becomes 

w(x,y,t)=A  ^/j(0)w(x,y,t)+J"^w(xi— VT,x2,r)/i'(t— r)  drj 

=A^0)w(x,y,t)+j^w(x+vt-VT,y,r)/i'(t-r)  drj. 

Application  of  the  Fourier  transform,  defined  by  f(p,y,t)=l  e1Pxf(x,y ,t )  dx, 

-*—ao 

to  equation  (2.5)  results  in 

t 

p[lf+ivp]  w(P>y-t)=  [(-ip)2+^2]  [/<0)w(p,y,t)+  ^w(p,y,T)e-ivp(t_rV'(t-r)  dr  . 

Subsequent  application  of  the  Laplace  transform  g(p,y,s)=  g(p,y,7')e-sr  dr 

-*0 

to  (2.6)  yields  the  equation 

p(s+ivp)2w(p,y,s)=  [|^f-P2]  [^s+ivp)w(p,y,s)j 

in  which  £(s)  is  the  Carson  transform  of  the  shear  modulus  given  by 

p(s)=s/i(s)=p(0)+ p  e-BTd (i. 

Jo 

Equation  (2.7)  can  be  rewritten  as  -j^w—  jp^+- — — (s+ivp)2|w=0 

L  £(s+ivp)  J 

which  has  the  solution  w(p,y,s)=A(p,s)e'_l^s’p^^  ^ ,  where 

/9(s,p)=  [p2H - £ — (s+ivp)2]  must  be  chosen  so  that  Re  /?>0. 

*-  £(s+ivp)  J 

In  a  similar  manner,  the  Fourier  and  Laplace  transforms  may  be  applied  to 
the  constitutive  equation  (2.4)  to  produce  <723(p,y,s)=$s+ivp)^w(p,y,s).  If  one 
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defines  f^p)  and  f  (p)  by  f"*”(p)=  f™ e'^xf(x)  dx  and  f  (p)=  e'Pxf(x)  dx,  then 

Jq  '’-00 


’0 

the  boundary  condition  (2.3)  transforms  to 

^2t(p,0,s)+<T2s(p)0)s)=/3(s+iFp)^w(p)0)s)=-/a(s,p)/i(s+ivp)w(p)0(s). 

From  (2.3)  it  can  be  seen  that  <r23(p,0,s)=LeaeAe(aep,s)=g(p,s)  and 

w(p,0,s)=w— (p,0,s).  It  is  assumed  a  priori  (and  is  easily  verified  a  posteriori)  that 

w~(p,0,s)  and  tf^tXpAs)  have  analytic  extensions  w~(z,0,s)  and  <72"£(z,0,s)  for 
lm(z)<0  and  lm(z)>0,  respectively,  which  vanish  as  |z|-»®.  Thus  the  transformed 
boundary  condition  (2.8)  can  be  recast  as  the  Riemann— Hilbert  problem:  find 
F"^(z)  analytic  for  lm(z)>0  and  F~(z)  analytic  for  lm(z)<0  such  that 
1  im  F"^(z)=  1  im  F— (z)=0  and  on  lm(z)=0, 

Im(  z  )-*+oo  Im(  z  )  -4—00 

F+(p)=T(p)F“(p)-g(p,s)  for  pe(— a,x) 

where  F+(z)=<72"3(z>0,s)>  F— (z)=w-(z,0,s),  and 

~  .  T  2  o  2l^2  •  . 

T(p)=— £(s-Hvp)  p  + — - (s+ivp)  .  For  convenience,  explicit  reference  to 

*-  p(s+ivp)  -I 

the  s  dependence  of  T(p)  and  F*(z)  is  being  suppressed.  It  suffices  to  note  that  all 
of  the  following  analysis  is  valid  for  any  complex  s  with  Re(s)>0. 

It  is  well  known  that  the  solution  of  (2.9)  is 


£  (T) 
where  X  (z)  solves  the  homogeneous  Riemann— Hilbert  problem 

X+(p)=T(p)X_(p)  for  pe(^»,a>). 

To  solve  (2.11),  it  is  convenient  to  factor  T(p)  into  the  product 


T(p)=TI(p)T2(p)T3(p)  in  which  T,(p)=-£(s+ivp),  T2(p)=  [{r^j2]  ' \  and 


(2.8) 


(2.9) 


(2.10) 


(2.11) 
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T*p) 


=  [-£- 
Hp~ 


-El. 


?  J3(s+i  vp> 


.  X  (z)  may  now  be  constructed  as  the  product 


X*(z)bX*(z)X2*(z)X3*(z)  with  each  Xj*(z)  satisfying  the  Riemann— Hilbert 
problem  Xi+(p)=Tj(p)Xi_(p). 

What  will  ultimately  be  required  is  to  solve  the  homogeneous 
Riemann— Hilbert  problem  (2.11)  for  each  fixed  s  on  a  Bromwich  path  with  Re(s)>0. 
This  will  be  accomplished  by  first  assuming  s  real  and  positive  and  then  invoking  an 
analytic  continuation  argument.  Additionally,  for  the  subsequent  analysis  the  shear 
modulus  will  be  assumed  to  be  positive,  continuously  differentiable,  non— increasing, 
convex  and  such  that  /j(oc)=lim  /i(t)>0.  Convexity  is  sufficient  but  certainly  not 

necessary  to  insure  the  validity  of  the  following  calculations  and  though 
theoretically  overly  restrictive,  it  holds  for  most  of  the  customary  models  such  as  a 
standard  linear  solid  or  a  power— law  material.  Moreover,  it  is  worth  noting  that  no 
explicit  time  decay  rate  for  the  shear  modulus  needs  to  be  specified  for  the  results  to 
be  valid.  From  the  fact  that  ^(t)=0  for  t<0,  it  easily  follows  that  [$s+ivz)]  *  is 
analytic  for  lm(z)<0.  Therefore  one  may  choose 
X[+(z)=l  and  X|~(z)=-[/I(s+ivp)]~1. 

Since  the  product  T2(p)T3(p)=/3(s,p),  the  branches  of  T2(p)  and  T^p)  must  be 
chosen  so  that  their  product  satisfies  the  requirement  that  Re  $s,p)>0.  This 
condition  can  be  met  by  choosing  the  branch  of  2  ^2  with  branch  cut  along  the 
negative  real  axis  for  both  T3(p)  and  T3(p).  (See  Fig.  2.1.)  Therefore  T3(p)  can  be 
expressed  as  T2(p)=sgn(p)(p— ^)  and  X2*(z)  may  be  chosen  to  be 
X2"*’(z)=w+(z)  and  X2~ (z)=uT(z)(z—  ^-)~* 

in  which  J^(z)=z^7  with  branch  cut  along  the  negative  imaginary  axis  and 
u/~(z)=z'^2  with  branch  cut  along  the  positive  imaginary  axis.  Finally,  one  may 


(2.12) 


(2.13) 


1 
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construct  X3*(z)  by 

X)±(J>=Mp(r1(z»  where  r4(2)=jU"  lo6gj(T)J  JT. 

J  “OC 

To  find  a  closed  form  expression  for  Xs^z),  it  is  necessary  to  determine  the 
mapping  properties  of  T3(z)  for  z  in  the  half— plane  lm(z)<0.  If  one  first  considers 
/I(s4-ivz)=^(s— vq+ivp)  on  the  horizontal  lines  z=p+iq,  q<0,  pe(— 00,00),  then  it 
follows  easily  from  the  stated  assumptions  on  ^j(t)  that 

(i)  #0)=/i(oo)</i(s-vq)<R,e  $s-vq+ivp)</i(0)=/!t(s-vq±ioo); 

(ii)  Im  £(s— vq+ivp)=— Im  £(s— vq— ivp); 

(iii)  arg  £(s-vq+ivp){>[J  ^  ; 

(iv)  1  im  /I(s— vq+ivp)=/i(0)  with  p(s—' vq)  converging  monotonically  to  ji(0). 
q->  -  00 

Therefore  £(s+ivz)  maps  these  lines  to  the  curves  shown  in  Fig.  2.2a. 

The  linear  fractional  transformation  $(z )~zl-i(s/v)  maPs  ^nes  z=p-t-iq, 
q<0,  pe(-oo,«)  to  the  circles  shown  in  Fig.  2.2b  and  S(p+iq)  exhibits  the  following 
properties: 

(v)  0=S(0)<S(0+iq)<Re  S(p+iq)<S(±oo+iq)=l; 

(vi)  Im  S(p4-iq)=— Im  S(— i3-f-iq); 

<’“)  «g  S(P+i<l)(<o  j£jj; 

(viii)  1  im  S(p+iq)=l  with  S(iq)  converging  monotonically  to  1. 
q-t  -  00 

If  the  elastic  shear  wave  speeds  corresponding  to  infinite  and  zero  time  are 
defined  by  c^=p(m)fp  and  c2=/i(0 )jp  then  from  (i)-(viii),  it  follows  that  (T^z)]2 
has  the  following  properties: 

(ix)  Im  T32(p+iq)=-Im  T32(-p+iq); 

(x)  arg  T32(p+iq)  P>§  ; 

(xi)  lim  T32(p-Hq)=l— (v/c)2; 
q-t  -  oo 


(2.14) 
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(xii)  -<v/c*)'fc<=^-=T32(0)<T32(iq)<T32(-4oo)=l-(v/c)‘'  where  T3“(iq)  is 

«s) 

monotonically  increasing  to  1 — (v/c)^  as  q-*-®. 

2 

Thus  it  can  be  seen  that  T3  (z)  maps  the  horizontal  lines  z=p+iq,  q<0,  pe(— ®,oo)  to 

>2 

the  curves  shown  in  Fig.  2.2c.  Furthermore,  it  can  be  seen  from  (xii)  that  T3(z) 

has  a  unique  root  z*=iq*  in  the  half— plane  lm(z)<0  for  any  positive  real  value  of  s. 

If  the  branch  cut  for  the  square  root  defining  Ts(z)  is  chosen  along  the 

negative  real  axis  then  T3(z)  is  analytic  for  lm(z)<0  except  for  the  branch  cut  on  a 

segment  of  the  negative  imaginary  axis  across  which  it  has  a  jump  discontinuity 

.2  _2  ,»/2 


given  by  T3(±0+iq)=±i 


L 

(q-|r  /J(s-vq)l 


for  q*<q<0.  It  follows  that  log  T3(z) 


is  analytic  for  lm(z)<0  away  from  this  line  segment  and  has  the  jump  discontinuity 
log(T3(+0+iq))— log(T3(— 0+iq))=ir  across  z=iq,  q*<q<0.  One  can  then  evaluate 


dr  as  in  Walton(1982)  and  conclude  that 

r1(2)=(l/2)[log(q,i-«Hog(-i)+log(Tj(®))  +  t 

Therefore  (2.14)  reduces  to 
Xj+(z)={W+(z-iq.)[T3(.)l'/2}/w+(2), 

From  equations  (2.12) — (2.15)  it  then  follows  that 
X+(z)=a,+(z-iq,)[T,(x)]'^, 

X— (z)=— {w~(z-iq#)[Ts(<»)]l/!}/[JZ(s+ivzXz-i(s/v))Tj(z)]. 

Finally,  for  a  specific  load  <r2s(*it)=LeAe(x/ae,t)  one  can  determine 

F"*’(z)=<r2^(z,0,s),  F— (z)=w~ (z,0,s)  from  equations  (2.10)  and  (2.16).  The  Laplace 
transform  of  the  stress  intensity  factor  (SIF)  K(s)  can  now  be  calculated  as  in 
Walton(1982).  In  particular,  it  is  straightforward  to  show  that 


(2.15) 


(2.16) 
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<723(x,s)~j^x  ^2  as  x-*0+,  where  K(s)=-[T3(ia)]^2  dr 

and  o'_(r)s)=LeaeAe(aer,s).  Again,  as  in  Walton(1982),  (2.17)  may  be  simplified 

to  K(s)=^-|  L  >/a  A  (x,s)jx|l^2e_^x^*ae^  dx,  in  which  it  should  be  emphasized 
v  J  —00 

that  q*  is  a  function  of  s  determined  implicitly  by  Tj(iq*)=0. 

Equations  (2.10)  and  (2.16)  can  be  shown  to  hold  for  each  fixed  complex  s 
such  that  the  Re(s)>0.  Therefore  these  expressions  (with  slight  modification)  are 
valid  on  a  Bromwich  path  and  numerical  Laplace  inversion  may  be  performed. 
These  details  are  found  in  Herrmann  and  Walton(19SS). 


L 


(2.17) 
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3.  Calculation  of  the  Energy  Release  Rate. 

The  Energy  Release  Rate  (ERR)  will  now  be  calculated  based  upon  the 

assumption  that  a  Barenblatt  type  failure  zone  exists  at  the  crack— tip.  Specifically, 

it  is  assumed  that  two  loads  are  acting  on  the  crack— faces:  the  applied  (external) 

tractions  denoted  <rT(x,t)=L  A  (|~,t)  and  the  cohesive  (failure)  stresses 

e 

<r^(x,t)=— L£A^w(x,t),^-,t)  acting  in  a  failure  zone  of  length  a^  immediately  behind 

the  crack— tip.  Two  essential  features  of  the  Barenblatt  model  are  that  a^<<ae  and 

that  Ke+K|=0  where  Kp  and  are  the  SIF'a  corresponding  to  and  crj, 

respectively.  The  effect  of  the  failure  zone  is  to  cancel  the  singular  stresses  ahead  of 

the  crack-tip  and  thereby  produce  a  cusp  shaped  crack  profile  behind  the  tip.  The 

resulting  mathematical  problem  is:  given  <r~(x,t),  find  the  "response"  stress  in  the 

failure  zone  and  crack  face  displacement  w(x,t)  that  cancel  the  stress  singularity  due 

to  *r~(x,t)  while  maintaining  a  constant  crack  speed.  The  goal  is  then  to  compute 

the  time  evolution  of  the  energy  release  rate  as  described  below. 

The  ERR,  G(t)  (defined  to  be  the  energy  flux  into  the  crack  tip  per  unit 

lfvt  -  .  . 

crack  advance)  is  given  by  G(t)=-  (t^xj— vt,t)u3(xi— vt,0,t)dxi,  which  in  the 

Jvt-a{ 

moving  coordinates  becomes 


G(0=4j  ^-t )  [|r -^]w(x,0,t)dx. 


Deriving  a  closed  form  expression  for  (3.1)  valid  for  arbitrary  applied  tractions 

a~ (x,t)  and  failure  zone  constitutive  law 

aj(x,t)=-L{A((w(x,t),|-,t) 

requires  solving  an  exceedingly  complicated  non-linear  boundary  value  problem. 
The  constitutive  law  (3.2)  models  the  failure  zone  response  as  that  of  a  nonlinear 
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elastic  spring,  accounting  for  inhomogeneity  and  aging.  In  general  the  failure  zone 
length,  aj,  is  dynamically  changing  and  determined  by  the  equation  Ke+K^=0.  In 
the  absence  of  any  generally  accepted,  physically  motivated  form  for  A^(  *,*,•). 
various  ad  hoc,  artificial  models  have  been  introduced  and  studied  in  the  literature. 

In  Schapery  (1975),  it  is  argued  that  taking  ^—(x,t)  constant  for  — a^<x<0 
with  aj=a£(v,t,<re  (x,t))  is  a  reasonable  approximation  to  reality,  at  least  provided 
a|  is  very  much  smaller  than  all  significant  physical  length  scales  such  as  crack 
length  and  the  distance  from  the  crack  tip  to  the  boundary.  More  recently, 
Goleniewski(1988)  studied  this  same  failure  zone  model  for  the  corresponding 
dynamic,  steady-state  problem  for  a  Maxwell  fluid  and  piecewise  constant  applied 
tractions  <r  “ (x,t).  However,  the  methods  of  Goleniewski(1988)  do  not  seem  to 
generalize  to  the  transient  problem,  even  for  a  Maxwell  fluid  Due  to  the 
complicated  non— linear  coupling  between  w(x,t)  and  aj,  deriving  a  convenient 
expression  for  G(t)  from  which  qualitative  properties  can  be  inferred  when  cr^  is 
assumed  constant  seems  unlikely. 

In  contrast,  for  the  dynamic  steady-state  problem  it  was  shown  in 
Walton(1987a)  that  a  simple  closed  form  expression  for  G  is  obtained  for  the  special 
class  of  loads 

<re(x)=Leexp(x/ae)  and  <rj(x)=— L^exp(x/a^)  for  -oo<x<0. 

It  was  argued  there  that  for  a£<<ag,  the  fact  that  0j(x,t)  does  not  have  compact 
support  should  have  a  relatively  minor  effect  on  the  results  provided  the  essential 
requirements  for  the  Barenblatt  model  are  still  satisfied:  a^<<ae  and  Ke+Kj=0. 
Furthermore,  the  two  cases  (1)  constant  with  aj=a^(v)  and  (2)  a^  constant  with 
Lf=L|(v)  were  compared  quantitatively.  Tt  was  found  that,  except  for  very  high 
crack  speeds,  the  two  cases  produce  nearly  identical  G  vs  v  curves.  In  light  of  these 


(3-3) 
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considerations  it  is  likely  that  valuable  insight  into  the  combined  influence  that 
crack  speed,  inertia  and  viscoelastic  properties  have  upon  the  rate  of  convergence  to 
steady— state  of  the  transient  G(t)  can  be  obtained  from  generalizing  the  analysis  of 
Walton(  1987a)  to  the  transient  problem  addressed  here,  the  unphysical  nature  of 
the  failure  zone  constitutive  law  (3.3)  notwithstanding. 

Henceforth,  the  forms  (3.3)  will  be  assumed  for  o g  and  Moreover,  time 
dependent  tractions  a~ (x,t)  will  be  permitted  by  taking 

<r“(x,t)=Lele(t)exp(x/ae),  -oo<x<0,  (3.4) 

where  1  (t)  is  a  dimensionless  function  of  time.  As  with  the  steadv-state  analvsis  in 

Walton(1987a),  consideration  of  the  two  models,  (1)  Lj  constant  with  a£=aj(v,t) 

and  (2)  a^  constant  with  Lj=Lj(v,t),  arises  naturally.  While  case  (1)  is  the  more 

physically  compelling,  it  is  mathematically  much  more  complicated  owing  to  the 

nonlinear  manner  in  which  the  initially  unknown  function  a^(v,t)  occurs  in  the 

problem.  In  contrast,  case  (2)  is  clearly  unphysical  in  a  dynamic  analysis  (e.g.  it 

suggests  that  information  travels  with  an  infinite  speed  of  propagation  in  the  failure 

zone)  but,  as  shown  below,  admits  an  elegant  closed  form  expression  for  G(t)  with 

the  aid  of  which  qualitative  and  quantitative  properties  can  be  easily  studied.  In 

light  of  the  steady  state  results,  it  is  likely  that,  except  for  very  high  crack  speeds, 

the  two  cases  should  exhibit  similar  long  time  asymptotic  behavior  whenever 

a|<<ag.  The  analysis  of  case  (2)  is  presented  here  and  that  of  case  (1)  is  the 

subject  of  a  future  paper. 

With  <re_(x,t)  given  by  (3.4)  and  <rf-(x,t)  by 

<rj(x,t)=Lflj<t)exp(x/a{),  -«<x<0,  (3.5) 

the  appropriate  definition  for  G(t)  becomes 

G(0=4j°  <^x,t)[|jr-v|^]w(x,0,t)dx. 


(3.6) 
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In  (3.5)  the  explicit  dependence  of  l^(t)  upon  v  has  been  temporarily  suppressed.  As 
shown  in  Walton  (1987a),  it  is  straight  forward  to  extend  the  analysis  to  treat  more 

general  loads  of  the  form  cr~~(x, t)s=L(t)  |  erx^a  dh(r)  where  h(r)  is  any  signed 

J0 

measure  for  which  the  integral  makes  sense.  However,  for  the  sake  of  brevity  that 
development  is  not  included  here, 

Incorporating  (3.5)  into  (3.6)  there  results 


G(t)=I(t)+W(t)  with 


SxlAf)d 


ri£w(x,0,t)  dx  and 


''-00 

W(t)=Lfl{(t)j<)  e(x/aP|rw(x,0,t)dx. 


I(t)  and  W(t)  in  themselves  do  not  have  any  direct  physical  significance  but  the 

decomposition  (3.7)  does  provide  some  insight  into  the  behavior  of  G(t).  In 

particular,  it  is  shown  below  that  lim  G(t)=lim  I(t)  and  lim  G(t)=lim  W(t) 

t-*0  t-*0  t-*oo  t-*oo 

whereas  lim  W(t)=lim  l(t)=0-  Thus  the  short  time  behavior  of  G(t)  is  governed 
t-*0  t—*  QC 

primarily  by  I(t)  and  the  long  time  behavior  by  W(t).  Moreover,  in  anticipation  of 
comparison  with  future  results  from  the  analysis  of  case  (1)  l^(t)=l  and  a^=a^(v,t), 
both  G(t)  and  W(t)  are  studied  in  detail  here. 

Motivation  for  the  separate  consideration  of  W(t)  derives  from  the  following 
observations.  In  Schapery(1975),  it  is  argued  that  it  is  physically  more  compelling 
to  build  a  viscoelastic  fracture  criterion  upon  the  work  done  to  the  trailing  ligament 
in  the  failure  zone,  T(t),  rather  than  upon  G(t),  the  total  energy  flux  to  the  failure 
zone  per  unit  crack  advance.  Of  course  this  concept  only  has  meaning  for  a  failure 
zone  with  finite  length.  For  any  failure  zone  constitutive  model,  one  may  write 


(3.7) 

(3.8) 

(3.9) 


20 


I 


G(t)  as  in  (3.7)  with  W(t)=j  ^  (x,t  )|_w(X)0,t)dx.  For  the  particular 

constitutive  model  (x,t)=— L^,  -a^<x<0  with  a^=a^(v,t),  studied  by  Schapery, 

r(t)=W(t).  It  then  seems  reasonable  to  study  W(t)  as  a  suitable  generalization  of 
r(t)  to  the  constitutive  model  (3.5)  for  case  (1)  lj(t)=l,  a£=a^v,t).  Expecting  case 
(1)  and  case  (2),  a^  constant,  lj=lj(v,t),  to  exhibit  similar  behavior,  at  least  as  t-*®,  a 
separate  analysis  of  W(t)  given  by  (3.9)  is  included  here. 

I(t)  and  W(t)  can  be  constructed  by  separately  evaluating  the  integral 
expressions  in  (3.8)  and  (3.9)  and  multiplying  the  results  by  l^(t),  which  is 
determined  from  the  Barenblatt  assumption  that  Kg+K£=0. 

Consider  first  the  integral 

g(t)4f  (3.10) 

If  one  notes  that  the  inverse  Fourier  transform  of  H(— x)e^x^af^  is  ^(p+i/a^)- ^  and 
applies  Parseval's  relation  to  (3.10)  and  then  applies  the  Laplace  transform  to  the 
resulting  expression,  the  Laplace  transform  g(s)  is  found  to  be 

[f  +*p]  w~(P'0>s)p4.f [£  •  Since  w~(P,0,s)=F~(p)  has  an  analytic 

extension  F— (z)  for  lm(z)<0,  the  integral  can  be  evaluated  by  residues,  whence  it 
follows  that 

+  (3ii) 

It  remains  to  evaluate  F~(— i/a^).  To  this  end,  one  begins  by  first  noting  that  (2.8) 
can  be  rewritten  as 

F~(p)=i(p)/T(p),  (3.12) 

where  a{p)=a~ (p)-t-ir+(p). 

For  the  Barenblatt  model  a~ (x)=<r~(x)-f  <^(x)  and  consequently 


(3.13) 
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a  (p)=^e(p)+^j(p)-  Also  cr+(p)  in  (3.13)  can  be  determined  by  application  of  the 
Plemelj  formula  to  (2.10)  thereby  obtaining 

<r"*’(p)=F"^’(p)=— ig(p)+ ^rX"*’(p)  P  -j-j.  This  can  be  rewritten  using 

*’-^30  X+(r)r-p 

^+(p)=^e(p)+^‘f(p) 

wherein  <r£(p)  and  <r"j(p)  are  given  by 


+/_x  lT-v_h,  1  v+,  X  ~*Jr)  dr 

a(p)=-^a(p)+^X  (p^-^——. 

#  .  T_  LJ  (S)  T_ 

From  (3.4)  it  can  be  seen  that  cj p)=^ ^  ijd~))  an<^  therefore  ae(z) 

meromorphic  with  pole  z=i/ap.  Since  X"^(z)  is  analytic  for  lm(z)>0,  (3.15)  can  be 
calculated  by  residues  with  the  result 


e  e  t  X+(i/ae)J 

Similarly,  <r^(p)w— <rr(p)fl  —  ^  (p) — 1  which  combined  with  (3.14)  and  (3.16) 
1  P  l  X+(i/a,)J 


yields 


J(p)=2Ll£l 


Lele(a  ) 


L  { i  (  s  )  1 

kp-(i/ae))X+(i/aJ'(p-(i/»f))X+(i/a{)J 
Equation  (3.17)  can  now  be  simplified  by  the  Barenblatt  hypothesis  Ke+K£=0.  For 
the  special  loadings  given  by  (3.4)  and  (3.5),  equation  (2.17)  yields 

K«(»)— tT*C-)],/*  e"/4  ^e'g(S)  and  K,(.)— (T*.)]'^  e*i''4  1^!!, 

X  (i/ae)  X  (i/af) 

which  under  the  Barenblatt  hypothesis  results  in  the  identity 

LeIe(s)  LfIf(s) 

X+(i/ae)  X+(i/af)‘ 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.1 8) 


22 


r 
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■ 
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Equations  (3.17),  (3.18),  (3.12)  and  the  observation  that  T(p)=X"*"(p)/X  (p) 
then  yield 


"*  if  2  (ae+a()  x+C/a  )  ' 


One  may  now  substitute  (3.18)  and  (3.19)  into  (3.11)  thereby  obtaining 


1W- 


VeWfs  i  X 

l'  «fl  Ue+a()  x+(i/a  }  1 


l+afq* 

f/2  1  1 

0 

1  pv  * 

1'~ae9* 

'  ?(s+v/af)| 

"  1  '  0 

(l+sa{/v)*  p(s+v/aj) 

which  when  combined  with  the  result  (2.16)  derived  previously  for  X*(z)  produces 
~LeMs^ae-,af)  _ l1+a tf*!'2  1  I  i  m2  \~'h 

In  order  to  present  the  results  in  a  nondimensional  form,  it  is  necessary  to 

introduce  certain  parameters.  First,  a  nondimensional  shear  modulus  is  defined  by 

p(t)=p  m(t/r)  where  p  =lim  p(t)  and  thus  iim  m(t)=I.  Also,  the  nondimensional 

t— » 00  t— 'CD 

parameters  7,  e,  a,  and  /3  are  defined  by  7=v/c*,  a=c*r/ae,  and  /?=— apq* 

where  0<7<c/c*,  c<<l,  a>0,  and  #>0. 

The  Laplace  transform  g(s)  can  now  be  rewritten  in  nondimensional  form  as 

— L  a  I  ( s L/f  ri 
e  e  e1-  /v  (I— e 


g(s)= 


27^ 


•(ra} 


i-<a 

V*  1  1 

(a/d2  1  I 

|-*/» 

m(rs+a7/e) 

(TS+a7/e)^  m(rs+a7/f)| 

Also  note  that  /fc=$(s,7)  is  the  unique  root  of  the  equation 
m(  7S+  <*7/?)= 7^(  1 +( Tsfay/3)  )2. 


Attention  will  now  be  turned  to  l^(t).  It  follows  from  (2.16)  and  (3.18)  that 

V8)=ys)r^|4T/?  • 

Finally,  utilizing  (3.7)— (3.10),  (3.20),  and  (3.21)  one  may  write  the  ERR  as 
L  2a 

G(t)=-L{l{(t)g(t)=:-^ril(t)gI(t), 


(3.19) 


(3.20) 


(3.21) 


(3.22) 
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where  the  Laplace  transforms  of  l(t)  and  g(t)  are  given  by 


gl(s)=- 


l+*/?l1/2 

~t+P\  ’ 

U-<)|1-<||,/2  1  j 

(q/02  1 

U+41  1+3  £(t s+07/e) 

2 

(rs+ati/e)  m(rs-fa7/e) 

'/  2 


Unfortunately,  it  is  not  possible  in  general  to  invert  (3.23)  and  (3.24) 
analytically.  Therefore  to  obtain  a  better  understanding  of  G(t)  either  an 
asymptotic  analysis  based  on  (3.22)— (3.24)  may  be  performed  which  is  valid  for  any 
shear  modulus  //( t)  or  numerical  inversion  of  the  Laplace  transforms  (3.23)  and 


(3.24)  may  be  obtained  for  a  specific  fj( t).  The  asymptotic  analysis  and  the 
consideration  of  special  cases  are  addressed  in  the  next  section.  Numerical  inversion 


of  (3.23)  and  (3.24)  for  the  standard  linear  solid  and  power— law  models  can  be 
found  in  Herrmann  and  Walton(1988). 


(3.23) 

(3.24) 


I  1 

1 
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4.  Special  Cases  and  Asymptotic  Solutions. 

Asymptotic  expansions  for  the  ERE,  G(t),  as  t-0  and  as  t-*x  can  now  be 
constructed  from  (3.22)— (3.24).  More  specifically ,  asymptotic  expansions  for  l(t) 
and  gift)  can  be  constructed  separately  and  multiplied  together.  Moreover,  it  can 

be  seen  that  l(t)=Ie(t)*l|(t)  where  Ii(s)=  •  Unless  indicated  to  the 

contrary,  henceforth  attention  will  be  limited  to  the  special  case  le(t)=2,  i.e.,  the 
external  tractions  driving  the  crack  will  be  assumed  to  be  time  independent.  In  that 
case  one  clearly  sees  that  I ( s)=s  ^Ii(s).  More  general  time  dependence  for  l^f  t )  is 
easily  incorporated  into  l(t)  by  a  simple  convolution  of  l^ft)  and  lift ). 

For  short  time  approximation,  it  is  necessary  to  determine  the  asymptotic 
behavior  of  I (s)  and  gi(s)  as  s-*oo.  To  this  end,  it  is  required  to  study  the  behavior  of 
rn(-ns+a7/e)  and  tfas  s-oo.  From  the  definition  of  the  Carson  transform  infs),  it  is 

o 

easily  seen  that  m(x)=m(0)=^^-=  —  .  Furthermore,  since  0=/9( s,j)  satisfies 

^oc 

the  equation 

m(7s+a7^)=72(l-f(rs/o'7/3))2  (4.1) 

it  can  be  shown  (see  appendix)  that 

+  °<s>  M  '4  2> 

Therefore  if  ifs)  and  gi(s)  are  written  as  asymptotic  series  in  powers  of  it 
is  found  that  as  s-*x 

I(s)=i/e[|  +  ^~^$(c/c*)-t]— < o  +  o(s-2)l,  (4.3) 

L  rs  J 

g«W-i 

Hence,  if  one  assumes  that  !(t)  and  gift)  have  Maclaurin  expansions  in  a 
neighborhood  of  t=0  then  from  (4.3)  and  standard  asymptotic  results  for  the 
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Laplace  transform, 


l(t)=V?[l  +  ^7^§(c/c*)-7](t/T)  +  oft)], 

Jm(V)  1  "  +  +  o(t)  35  w0' 


"  pT^0^-7]  +  +  o(t)  35  w0- 

From  (4.4)  and  (4.5)  it  follows  after  algebraic  simplification  that 
L  2a 

v f  [l  -aWT)  +  °m]  15  ^°> where  a=ie(c-v)  +  TH^' 

It  should  be  remembered  that  (4.4)  and  (4.5)  (and  hence  (4.6)  also)  are  not 
uniformly  valid  expansions  as  e-+0.  The  case  e=0  is  a  singular  limit  treated 
separately  later.  It  should  also  be  remarked  that  near  t=0.  G(t)  is  governed  by  the 
glassy  properties  except  for  the  inclusion  of  the  term  m'(0)  which  incorporates  the 
effect  of  the  initial  rate  of  stress  relaxation.  Since  m(0)>l  and  m'(0)<0,  one  sees 
that  the  sign  of  G'(0)  depends  upon  the  crack  speed  and  material  properties  through 
the  relation 

G'(0)>0  |G'(0)<o]  if  and  only  if  ae^g|i->[c-v]  {aJ^|jjU[c-v]}. 

In  particular  for  fast  enough  crack  speeds,  G'(0)  will  be  positive,  i.e.  G(t)  will 
initially  increase  with  time.  However,  for  any  given  crack  speed  if  the  combination 
^  '5  sma^  enou&l)'  then  G(t)  will  initially  decrease  with  time. 

To  determine  long-time  asymptotic  solutions  as  t-*oo,  it  is  necessary  to  find 
L  2ae 

first  1  i  m  G(t  )=-£ — -  1  i  m  si (s)  1  i m  sgi(s).  The  expression  m(  rs+  07/ e )  has  the 
t-*  00  ^00  s-*0  s-*0 

limit  mfcry/e)  as  s-*0  and  $  as  shown  in  the  appendix,  has  different  asymptotic 
limits  as  s-*0  depending  on  whether  (i)  0<v<c*,  (ii)  v=c*,  or  (iii)  c*<v<c,  namely 

+  °(s)  *or  ®<v<c*i 

rra'(r)  dr]  ^s^2  +  o(s'/2)  for  v=-c*, 

Q2L  J0  J 

P=PQ  +  o(l)  for  c*<v<c,  where  /3q  is  defined  implicitly  through  the  equation 
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7^—  ji--" J  =  J^e  a^orm'(t)  dr.  To  avoid  separate  cases  in  displaying  subsequent 

formulas  it  is  convenient  to  define  B  =0  for  0<v<c*. 

o  ~ 

Again  from  standard  asymptotic  arguments  for  the  Laplace  transform 

Le2\n-,)  1  T2  -‘h 

which  is  precisely  the  steady-state  solution  found  in  W alt on(  1987a), 

It  is  interesting  to  consider  separately  the  asymptotic  behavior  of  L^t)  for  a 
general  time  dependence  LeI^(t )  in  cr“(x,t ).  One  notes  that  as  t-*0, 

Lflj(0)=lim  sI/j!f(s)=Lele(0).  Thus  initially  L^t)  and  Lele(0  are  equal  for  any 

crack  speed  and  f>0.  In  steady— 6tate,  however. 


L£l^oo)=lim 
!  s-»0 


d/0+/?o 


xh 


.  For  0<-v<l,  ^  =0,  which  implies  that 


L^l^f  «)=( l/>/7)Lele(Qo)  which  is  the  same  as  found  for  quasi— static  crack  propagation. 

For  l<7<c/c*,  /?o  is  a  monotone  increasing  function  of  7  with  lim  /?0(  7)=*- 

T4  c  / c* 

Lflf(oo) 

Thus  ^  is  a  monotone  decreasing  function  of  7  with  1  i  m  L^oo)=Le(oc).  For 
e  e'”'  7-*c/c* 

the  special  case  of  L  I  (t)  constant,  one  obviously  uses  L  1  (0)=L  1  (oo)=L  in  the 
€  e  c  c  c  c  c 

above  formulas. 

It  is  instructive  to  consider  several  special  cases.  Setting  e=0  corresponds  to 
a  model  with  no  failure  zone.  Allowing  e-*0  in  (3.23)  and  (3.24)  yields 

Ks)=ie(s)A/T+?  and  gi(s)=Ie(s)^-gy  Jfirp-  [M*/c)2]  ^2-  It  then  follows  easily 
from  (3.18)  and  (3.22)  that  Gnj(t),  i.e.  G(t)  based  on  the  singular  stress  field 
without  inclusion  of  a  failure  zone,  becomes 
G^»=!K(t.7))2[^J|l-(v/c)2)'',J 


(47) 


(4.8) 
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Thus,  Gn|(t)  is  simply  the  product  of  the  square  of  the  SIF,  K(t,7)  ,  with  a  time 
constant  function  of  crack  speed  and  the  glassy  material  properties. 

If  one  assumes  that  K(t,7j  has  an  asymptotic  expansion  in  powers  of  t*^2  a 
t— *0,  then  it  can  be  shown  from  (2.16),  (3.18),  and  (4.2)  that 


K(t,7)=— LeVtaea/T)t(c/c*j 

produces 

L.^a. r  1  / 2 


l"r  +  o (Jt)  as  t-*0.  Substitution  into  (4.8)  then 


!(lc/ae)  +  °W  “  lH°- 

The  steady— state  limit  when  t-too  is  found  to  be 

It  should  be  remarked  that  the  right  hand  side  of  (4.10)  can  be  shown  to  vanish  as 


The  special  case  of  an  elastic  material  is  obtained  from  the  general  solution 
(3.22)-{3.24)  by  letting  r-*a>.  There  results 
l  2a 

G(l)=-T^|(t)gi(t), 

'  IH'/c  +  s*i/c)2r'^  wh«e 

The  quasi— static  limit  is  obtained  from  (3.23)  and  (3.24)  by  letting  c*-*oo 
(p-*0).  It  follows  that  0=0  and  hence  that 

T(i)=Te(s), 

In  particular,  one  sees  that  1^ f l|( t ) = Lele( t ) / -^e  and  thus  that  L^t)  is  crack  speed 
independent.  Moreover,  consideration  of  (4.12)  reveals  that  gi(t)  in  the  quasi— static 


(4.10) 


(4.11) 


(4.12) 
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case  has  a  Dirac  delta  function  singularity  at  t=0. 

Several  comments  on  these  asymptotic  results  are  in  order.  First,  it  should 
be  observed  that  whether  or  not  the  model  assumes  a  failure  zone  dramatically 
affects  both  qualitatively  and  quantitatively  the  behavior  of  the  ERR  as  a  function 
of  time,  crack  speed,  and  material  properties.  For  example,  from  (3.20)  it  is  easily 

shown  that  for  any  t,  if  e>0,  1  im  G(t)=®  whereas  if  6=0,  1  im  G  > 

7^0  r*  0  nl 


-Iv/aJe^ 


where  K(s,0)=  and  from  (4.1)  /3fs,0)=2-^[m(rs)]1^2.  However,  in  the 

steady— state  limit  (4.7),  G  approaches  a  finite  limit  as  the  crack  speed  vanishes, 
both  with  and  without  a  failure  zone.  Specifically ,  it  is  easily  seen  that  for  the 
steady-state  limit 

lim  G(oc)= — ry—  ~  for  e>0 

' rO  1 

L  2a  . 

— — -  e  -  for  6=0. 

2  M0 

Thus  G  becomes  infinite  as  7-O  except  for  6=0  (no  failure  zone)  or  under 
steady— state  conditions. 

The  reason  for  this  behavior  is  found  through  consideration  of  the  c rack  face 

particle  velocity  U3(xi,0,t)=-^-(x,0,t)— v^(x,0,t).  A  consequence  of  the  assumption 

that  there  is  an  initial  jump  discontinuity  in  the  applied  crack  face  tractions  is  that 

^j-(x,0,t)  does  not  vanish  as  v-*0.  Thus,  from  (3.6)  it  follows  that  1  im  G(t)=®.  In 

v-*  0 

2 

contrast,  when  6=0,  one  sees  from  (4.8)  that  G  is  merely  a  product  of  K(t)  and  a 
simple  function  of  crack  speed  and  glassy  material  properties  that  is  independent  of 
the  crack  face  particle  velocity  and  that  remains  bounded  as  v-*0.  Moreover,  in 


rO 


steady— state,  G  is  given  by  G=—  <r^  (x)^(x,0)  dx  and  thus  remains  bounded  as 


—00 
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v-*0. 


Other  differences  between  the  e=0  and  e>0  cases  are  evident  in  the  short 

time  behavior  of  G(t).  in  particular,  from  (4.6)  it  is  easily  seen  that  when  e>0, 
r  2 

L  a  , _ 

lim  G(t)=— ^  —  y+j  j  f>0  whereas,  for  e=0,  it  follows  from  (4.9)  that 

lim  G(t)=0. 

t-0 

The  singular  behavior  for  G(t,v)  as  v-*0  in  the  dynamic  solution  and  the  delta 
function  contribution  at  t=0  in  the  quasi-static  solution  are  both  due  to  the 
integral  I(t)  in  (3.7)  which  contains  the  relative  crack  face  particle  velocity 
■^•w(x,0,t).  As  shown  below,  W(t)  exhibits  none  of  the  singular  behavior  seen  in 

G(t). 

An  expression  for  W(t)  is  easily  constructed  from  simple  modifications  of 
(3.10),  (3.11),  and  (3.22).  Specifically,  one  sees  that 
L  2a 

W(t)=-r^Wg2(i) 


where  l(s)  is  given  by  (3.20)  and 


'/*  i 

o 

1  (  7+s  t  e  /  a) 

1+^ 

m(rs+a7/f) 

m(  rs+try/e) 

— '/2 


It  is  trivial  to  show  that  the  limit  of  W(t)  as  e-*0  (no  failure  zone)  is  Gnf(t)  given  in 
(4.8).  The  other  special  cases  considered  above  are  easily  constructed  for  W(t).  In 
particular,  for  elastic  material  l(s)  is  still  given  by  (4.11)  whereas  gi(s)  is  replaced 
by  g2(s)  given  by 

g2(s)=le(s)j^j|^3  ^  [Mv/c  +  saf/c)2]  '/2,  with  0=^,  as  before. 

The  quasi-static  limit  of  W(t),  W  (t),  (c*-*«)  yields  as  before  T (s)=  1  (s) 

os  c 

and  has  gi(s)  replaced  by  g2(s)=l  (s)fc*<  - - - - .  In  particular,  if  1  (0=1, 

^  r  '  m(  Ts+vr/aj) 


(4.13) 


(4.14) 
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i.e.  the  applied  tractions  are  time  independent,  then  I  (s)=-  and  it  follows  that 

6  S 

ga(t)=|y~j  m(TJ)  +  J  e- ^rvT^aPc'(r)  drj,  where  c(t)  is  the  non-dimensional 


compliance  corresponding  to  m(t)  defined  by  m(s)=l/c(s).  Thus, 

2 

_  ,  .  H—  f  1  _/rtTT/a.^  ...  .1 


L  2a  n 


£fc  +  .f 


It  should  be  noted  that  1  im  W  (t)=L  ^a  D(t/r),  where  D(t/r)  is  the  creep 

v-*0  qs  e  e 

compliance  corresponding  to  the  shear  modulus  /i(t/r). 

Short  time  asymptotic  expansions  for  W(t)  can  be  constructed  as  before  for 
G(t).  In  particular  from  (4.4)  and  (4.13)  it  follows  (assuming  1  (t)=l)  that 

L  2a 

w(t)=~2~ [x+t] r+°(t)as  t'4°+' 

A  similar  asymptotic  expansion  of  (4.14)  reveals  that  the  short  time  response  for 

W(t)  in  (4.16)  is  unchanged  for  an  elastic  material.  Obviously  though,  the 

asymptotic  series  for  elastic  and  viscoelastic  materials  differ  in  their  higher  order 

terms.  Moreover,  as  noted  previously  for  G(t),  the  asymptotic  expansions  for  l(t) 

and  g2(t)  are  not  uniformly  valid  as  e-*0.  Thus,  the  limit  of  (4.16)  as  e-0  does  not 

yield  the  small  t  expansion  (4.9)  for  Gn^(t).  It  should  also  be  noted  that  the  first 

term  in  the  asymptotic  expansion  (4.16)  is  independent  of  crack  speed.  Thus,  to 

leading  order  as  t-*0,  W(t)  for  a  moving  crack  is  the  same  as  that  for  a  stationary 

crack.  However,  the  higher  order  terms  are  crack  speed  dependent  and  influence 

W(t)  for  all  but  asymptotically  small  times  t. 

It  is  interesting  to  compare  the  dynamic  result  (4.16)  with  the  quasi-static 

result  (4.15).  In  particular,  the  dynamic  analysis  shows  that  lim  W(t)=0,  whereas 

t-*0 

L 

Le  ae  (1— el 

the  quasi-static  result  has  lim  W  (t)= — rj - fn  4.  Thus,  for  short  times  the 

qs  + 


(4.15) 


(4.16) 
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quasi— static  and  dynamic  analyses  yield  quite  different  results  for  all  crack  speeds. 
This  is  in  marked  contrast  to  the  steady— state  limits  which  show  good  agreement 
for  all  crack  speeds  up  to  nearly  (and  sometimes  beyond)  the  equilibrium  shear  wave 
speed.  The  short  time  asymptotic  analysis  reflects  the  obvious  fact  that  the  initial 
transient  response  for  the  quasi— static  case  should  differ  from  that  occurring  when 
the  effect  of  material  inertia  is  included. 

Attention  will  now  be  directed  toward  describing  the  manner  in  which  W(t) 
converges  to  the  steady— state  limit  as  t-»®.  It  is  useful  for  comparison  purposes  to 
consider  the  question  first  for  elastic  material.  Again  it  is  assumed  that  1  (t)=l,  or 
equivalently  I  (s)=s—*.  Thus, 

L  2a 

with  I(s)  and  g2(s)  given  by  (4.11)  and  (4.14),  respectively.  It  is  instructive  to 


consider  first  the  elastic  case  with  no  failure  zone.  In  that  case  l(s)=- 


1  1 


s  71+7? 


g2(s)=j^^l— (v/c)2]--1^2  with  0=  It  follows  easily  that  l(t)=Erf(7t/b) 
where  b=a e/(c— v)  and  hence  that  W(t)=-^~{l-(v/c)2]  ^2[Erf(7t/h)]2.  It  is 


easily  seen  that 


W(t)=W(oo)  +  0(e~(t/b))  as  t-®. 


Thus  W(t)  converges  exponentially  to  its  steady— state  limit  with  exponential  order 
e~^b^.  It  should  be  noted  that  b-*®  as  v-*c  which  implies  that  the  exponential 
decay  rate  vanishes  as  the  crack  speed  approaches  the  shear  wave  speed.  The  case 
of  an  elastic  material  with  a  failure  zone  is  slightly  more  complicated  but  can  be 


treated  in  a  similar  fashion.  One  has  I(s)=i-  and 

^s)=?[lT^a^^1+a)^15",/2tHv/c)2r‘/2  with  a=  ( .  A  simple 
calculation  then  yields  that  Wg^(t)  is  given  by  (4.17)  with 


(4-17) 


(4.26) 


1 
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U« wa„(d,t  +  B^jV(dIt)e^l'r  dr  and 

^‘'=bia!l+j]l1-('/c,2i  h\^'r  dr  where  -]=7T£  C2=2sS'  d|=27E' 

d2=5eF’  and  IQ(t)  is  the  modified  Bessel  function  of  the  first  kind  of  order  zero 
defined  by  I 

decay  exponentially  to  their  steady— state  limits  at  the  order  with 


b=ae/(c— v)  as  in  the  no  failure  zone  case.  Therefore  (4.18)  is  valid  also  for  an 
elastic  material  with  a  failure  zone. 


It  is  readily  apparent  that  for  a  viscoelastic  material  the  situation  is 

considerably  more  complicated  due  to  the  combined  influence  of  material  inertia  and 

viscoelastic  stress  relaxation  upon  convergence  to  steady  state.  A  general  property 

of  the  Laplace  transform  F(s)  of  a  function  f(t)  is  that  f(t)  decays  exponentially  in 

time,  say  f(t)~e~ 3  with  a>0,  if  and  only  if  F(s)  is  analytic  in  the  halfplane 

Re(s)>— a.  The  expressions  (3.23)  and  (3.24)  for  I(s)  and  gi(s)  were  derived  for  real 

positive  s.  Determining  the  largest  value  of  a,  &  ,  for  which  s I (s)  and  sgi(s)  are 

max 

analytic  in  the  halfplane  fte(s)>— a  is  a  difficult  task  that  clearly  depends  upon  the 
particular  details  of  the  transform  m{rs+ay/ e).  From  the  previously  stated 
properties  assumed  for  m(t)  it  follows  that  (fn(Ts+07/e)]— ^  is  analytic  for 

^  m(t)  is  a  power— law  in  t,  then  m(  Ts+a7/e)  has  s=~^  as  a  branch 
point  and  amax  can  be  no  larger  than  that.  As  a  second  example,  for  a  standard 
linear  solid  with  m(t)=l+ fje- *,  [m(  rs+ay/f)]"-^  is  analytic  for 
Re(  s  )>—((??+ 1  )_1 + <*yf  f]  /  t. 

The  heart  of  the  matter  lies  in  determining  the  analyticity  properties  of 


/3(s,7),  which  is  defined  implicitly  through  equation  (4.1).  Whether  or  not  s,7)  is 
analytic  in  some  halfplane  Re(s)>-a,  a>0,  depends  upon  the  particular  way  m(t) 
decays  to  its  equilibrium  value  as  t-*oo.  This  necessitates  a  case  by  case  analysis  for 
different  forms  of  m(t).  We  content  ourselves  here  with  illustrating  the  differences 
that  exist  between  materials  with  exponentially  decaying  modulus,  such  as  a 
standard  linear  solid,  and  those  for  which  m(t)  decays  as  a  power  of  t  to  its 
equilibrium  value,  such  as  a  simple  power— law  material  with  m(t)=l+rj(l+t)  n, 
n>0. 

An  important  observation  to  be  made  for  power— law  material  is  that  G(t) 

cannot  have  exponential  decay  to  its  steady-state  limit  when  7<1,  i.e.  v<c*.  The 

reason  for  this  is  that  /3(s,7)  is  not  analytic  at  s=0  as  can  be  seen  from  the  following 

argument.  It  was  remarked  earlier  that  /5(0,7)=0  whenever  7<1,  in  particular 

lim  /5(s,7)/s=-7-T^rT.  Moreover,  for  power— law  material  m(s)  has  s=0  as  a  branch 
s-*0  +  ^  ' 

point.  If  /5(s,7)  were  analytic  at  s=0  then  the  right  hand  side  of  equation  (4.1) 
would  also  be  analytic  there.  However,  the  left  hand  side  of  (4.1)  has  s=0  as  a 
branch  point.  This  contradiction  proves  the  claim. 

On  the  other  hand,  for  l<7<m(0)'^2  (i.e.  c*<v<c),  since  $0,7)>0,  equation 
(4.1)  defines  $3,7)  as  an  analytic  function  in  a  neighborhood  of  s=0.  In  particular, 
/5(s,7)  is  analytic  in  a  halfplane  Re(s)>-a,  a>0.  However,  finding  the  largest  such  a 
is  difficult.  For  elastic  material,  the  rate  of  exponential  convergence  of  G(t)  to  its 
steady-state  limit  corresponds  to  the  negative  read  value  of  s  for  which  7)= — 1 . 
For  viscoelastic  materials,  again  a  case  by  case  analysis  will  be  required  to 
determine  precisely  where  the  singularity  of  sl(s)  or  sgj(s)  with  largest  real  part  will 
occur. 

The  remaining  case  7=1  (i.e.  v=c*)  is  easily  handled.  Indeed,  for  a  general 
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material,  not  just  power— law  or  a  standard  linear  solid,  equation  (4.1)  admits  no 
solutions  for  s<0.  Thus  /3(s,l)  is  not  analytic  in  any  halfplane  Re(s)>--a,  a>0  and 
G(t)  cannot  converge  exponentially  to  its  steady-state  limit.  Furthermore,  it  is  not 
difficult  to  see  from  equation  (4.1)  that  the  left  endpoint  of  the  largest  s— interval 
containing  zero  for  which  an  admissible  solution  /9(s, -y)  exists  converges  to  zero  as  7 
approaches  1  from  above. 

These  observations  made  for  power— law  material  whenever  l<7<m(0)l/^2 
(c*<v<c)  are  equally  valid  for  a  standard  linear  solid.  A  departure  in  behavior 
occurs  for  0<7<1.  If  m(t)=l+^e— \  then  m(s)=^t|_~^^s  and  (4.1)  can  be 

J-'T  S 

rewritten  as 


where  z=T5+c»7/5(s,7).  When  s>0,  a  root  z  must  be  sought  that  is  greater  than  7s, 

whereas  z<rs  is  required  when  s<0  is  suitably  near  zero.  An  examination  of  the 

graphs  of  the  functions  on  either  side  of  equation  (4.19)  quickly  reveals  that 

admissible  solutions  exist  for  any  s>0  and  s<0  suitably  near  zero.  Moreover,  no 

solution  exists  for  7=1  and  the  left  hand  endpoint  of  the  largest  s  interval 

containing  zero  on  which  a  solution  exists  tends  to  zero  as  s-*0— .  One  concludes 

from  this  that  for  a  standard  linear  solid,  G(t)  converges  exponentially  to  its 

steady— state  limit  for  0<7<1  and  l<7<m(0),//2  but  not  for  7=1.  Moreover,  the 

rate  of  exponential  convergence,  i.e.  amax>  tends  to  zero  as  t-+1±.  It  is  also  easy  to 

see  from  (4.1)  that  a  must  varnish  as  T^m(O)1^2. 

'  '  max  ’ 


(4.19) 
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5.  Conclusions 

The  principal  contributions  of  this  paper  to  the  study  of  transient  mode  III 
crack  propagation  are  the  solutions  for  the  displacement,  stress,  and  stress  intensity 
factor  for  general  loadings  and  general  shear  moduli,  the  inclusion  of  a  failure  zone 
into  the  model  and  the  calculation  of  the  total  energy  flux  into  the  failure  zone, 
G(t).  It  was  then  observed  that  significant  qualitative  and  quantitative  differences 
exist  in  the  behavior  of  G(t)  as  a  function  of  time,  crack  speed,  and  material 
properties  between  a  model  incorporating  a  failure  zone  and  one  which  does  not. 

The  question  of  the  rate  of  convergence  of  G(t)  to  its  steady— state  limit  as 
t-*a  was  also  investigated.  It  was  observed  that  this  rate  of  convergence  depends  in 
a  complicated  way  upon  the  rate  of  stress  relaxation  and  crack  speed.  In  particular, 
for  a  standard  linear  solid  in  which  stresses  relax  exponentially  fast,  G(t)  converges 
exponentially  fast  for  all  crack  speeds  except  the  equilibrium  shear  wave  speed. 

Thus  for  crack  speeds  near  the  equilibrium  shear  wave  speed,  it  is  expected  that 
steady— state  conditions  would  set  in  more  slowly  than  for  crack  speeds  above  or 
below  it.  Also  the  exponential  rate  of  convergence  is  lost  at  the  glassy  shear  wave 
speed.  In  contrast,  for  power— law  material,  an  exponential  rate  of  convergence  of 
G(t)  does  not  occur  for  any  crack  speeds  less  than  or  equal  to  the  equilibrium  shear 
wave  speed  whereas  for  speeds  between  the  equilibrium  and  glassy  shear  wave 
speeds  G(t)  does  converge  to  steady  state  at  an  exponential  rate. 


V., 
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Appendix 

The  behavior  of  /?as  s-*®  is  determined  by  equation  (4.1).  If  one  lets 
h=/?/s>0  then  (4.1)  becomes 
m(  s[  r+ aTh] ) = 72(  1 +(  t/ 07b )  )2. 

The  left-hand  side  m(s[r+a7h])=m(0)+  f* e~^  T+Q^)srm'(r)  dr  has  the  limit 

Jo 

rc  i2  2 

m(0)=l— I  >7  ass-*®.  Therefore  (A. 1)  will  be  satisfied  only  if  h=h(B+o(l )  where 
h  >0.  As  s-*®,  (A.  1)  becomes  —I  =  72(l+(  7/07)1  ))2.  It  is  easily  seen  that 

00  |C  $  J  OP 

h.=51fc7cVP7) and  ,hus  ^olCc/^H-rl  +  °« as  ,H*- 

To  determine  the  behavior  of  /?  as  s-*0,  one  must  again  consider  equation 
(4.1).  Note  that  the  left-4iand  side  of  (4.1)  is 

m(  7s+a7/?)=m(0)+  f° e-^  ra+a^rm'(r)  dr>l. 

J0 

2  0 

For  0<v<c*,  the  limit  as  s-*0  of  the  right-hand  side  of  (4.1),  7  (1+ 75/07/?)  ,  will 
satisfy  (A. 2)  only  if  /3=/?|S+o(s)  as  s-*0.  Therefore  as  s-*0,  (4.1)  becomes 

l=7^(l  +  r/o7/?i)~. 

The  solution  of  (A. 3)  for  /?i  is  easily  seen  to  be 

For  c*<v<c,  the  right-hand  side  of  (4.1)  satisfies 
7^(  1+  t&!  07/?)2>  72>  1 . 

If  /M)  as  s-*0  then  it  can  be  seen  from  (A. 2)  that  the  limit  of  m(rs+a7/?)=l  as  s-*0 
which  contradicts  (A. 4).  Thus  /?=/?o-fo(l),  /3q> 0  as  s-*0  where  fiQ  satisfies 

^W-re^-W  dr- 

Jo 

For  v=c*,  it  can  be  shown  that  to  satisfy  (4.1),  /M)  and  s//M)  as  s-*0. 
Therefore  consider  /?s=/?is^+o(s^),  0<7<1,  as  s-*0.  If  one  substitutes  this  into  (4.1) 


( A.l) 


(A.2) 


(A.3) 


(A. 4) 
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and  simplifies  the  equation,  one  obtains 


rtO  L~m’Tf  -»  rOC 

— rm'(r)  dr  +  o(s^)=— —  +  c^s1-^)  as  s-*0,  assuming  that  rm'(r)  dr 

Jn  a7Pi  Jn 


■CD 

>0  ",y]  ^0 
exists.  If  this  integral  does  not  exist  then  one  must  do  the  asymptotics  for  the 

particular  shear  modulus  needed.  Since  the  coefficient  of  each  side  is  non-^ero, 
>7=1/2  and  j* rm'(r)  drj  .  In  summary,  4-  o(s)  for  0<v<c*, 

,^==^7[— rm'(r)  ^r]  s'^2  +  °(s^2)  f°r  v=c*i  d=,5o+o(l)  for  c*<v<c. 


I  1 

I  ft 


ACKNOWLEDGEMENT 


The  second  author  wishes  to  express  his  gratitude  to  Professor  J.R. 
Whiteman  and  the  Brunei  Institute  of  Computational  Mathematics,  Brunei 
University,  for  their  hospitality  and  support  during  the  sabbatical  year  in  which 
part  of  the  research  for  this  paper  was  conducted. 


39 


Achenbach,  J.D.  and  Chao,  C-C. 

REFERENCES 

1962 

J.  Mech.  Phys.  Solids  10,245 

Atkinson,  C. 

1979 

Arch.  Mech.  Stos.  31,829 

Atkinson,  C.  and  Coleman,  C.J. 

1977 

J.  Inst.  Maths  Applies  20,85 

Atkinson,  C.  and  List,  FL.D. 

1972 

Int.  J.  Engng  Sci.  10,309 

J.  Mech.  Phys.  Solids  27,431 

Atkinson,  C.  and  Popelar,  C. 

1979 

Baker,  B.R. 

1962 

J.  Appl  Mech  29,449 

Barenblatt,  G.I. 

1962 

Advances  in  ADDlied  Mechanics 

Goleniewski,  G. 

1988 

Vol  VII  Academic  Press.  55 

Int.  J.  Frac.  37,R39 

Herrmann,  J.M.  and  Walton,  J.R. 

1988 

Elastic-Plastic  Failure  Modelling  of 

Knauss,  W.G. 

1973 

Structures  With  Applications, 

Eds.  D.  Hui  and  T.J.  Kozik. 

ASME,  PVP-Vol  141,1 

Deformation  and  Fracture  of 

Popelar,  C.  and  Atkinson.  C. 

1979 

High  Polymers.  Eds.  H.H. 

Kausch  and  R.  Jaffe,  Plenum 

Pjggg  501 

J.  Mech.  Phys.  Solids  28,79 

Schapery,  R.A. 

1975 

Int.  J.  Frac.  11,141 

Schovanec,  L.  and  Walton,  J.R. 

1987a 

Acta  Mech.  67,61 

Schovanec,  L.  and  Walton,  J.R. 

1987b 

J.  Engr.  Frac.  Mech.  24,445 

Schovanec,  L.  and  Walton,  J.R. 

1987c 

submitted  for  publication 

Walton,  J.R. 

1982 

Quart.  Appl.  Math.,  April, 37 

Walton,  J.R. 

1983 

Workshop  on  Dynamic  Fracture. 

Walton,  J.R. 

1985 

Eds.  W.G.  Knauss,  K. 

Ravi— Chandar  and  A.J.  Rosakis. 
Cal.  Tech. ,197 

J.  Appl.  Mech.  52,853 

Walton,  J.R. 

1987a 

J.  Appl.  Mech.  54,635 

Walton,  J.R. 

1987b 

submitted  for  publication 

Willis,  J.R. 

1967 

J.  Mech.  Phys.  Solids  15,229 

I 


r- 


ft(s,o) 


